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Abstract 

We introduce a new effect in tfie neutrino oscillation phase which shows the 
neutrino— antineutrino oscillation is possible under gravity even if the rest masses of the corre- 
sponding eigenstates are same. This is due to CPT violation and possible to demonstrate if the 
neutrino mass eigenstates are expressed as a combination of neutrino and antineutrino eigenstates, 
as of the neutral kaon system, with the plausible breaking of lepton number conservation. For 
Majorana neutrinos, this oscillation is expected to affect significantly the inner edge of neutrino 
dominated accretion disks around a compact object by influencing the neutrino sphere which con- 
trols the accretion dynamics, and then the related type-II supernova evolution and the r-process 
nucleosynthesis. On the other hand, in early universe, in presence of various lepton number vi- 
olating processes, this oscillation, we argue, might lead to neutrino asymmetry which resulted 
baryogenesis from the B — L symmetry by electro-weak sphaleron processes. 
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I. INTRODUCTION 



It is well known that the discrepancy between various observed and expected number of 
neutrino flavors, e.g. i^e, and occurs because they transform from one type of flavor to 
another, namely flavor oscillation, which conserves lepton number. Therefore, if a detector 
is built and set to detect for one type of neutrinos, while that type of neutrinos are coming, 
some of them convert to another type which is unsuited to the detector and thus results 
a shortfall. This anomalies have been found in the cases of solar neutrino, atmospheric 
neutrino and the neutrino experiments at accelerators or reactors. This is also known that 
the neutrino flavor states consist of superposition of different mass eigenstates which evolve 
differently with time and thus flavor states transform from one to another. However, if CPT 
is violated, even for a particular flavor the neutrino and the anti-neutrino state may evolve 
differently and then the lepton number violating oscillation may take place. 

The study of neutrino has a lot of cosmological and astrophysical implications. For 
example, neutrinos are the important ingredient of the energy produced and removed from 
the center of Sun, they carry off the largest amount energy of an exploding star in supernova, 
neutrinos in certain accretion disks around a compact object are responsible for cooling 
process. Moreover, neutrinos may be considered as prime dark matter candidate. Although 
the mass of a neutrino is very less, because of their large number, density is very great 
and they would thus exert a very potent gravitational effect helping to determine the rate 
and the pattern of galaxy formation. Therefore, they may have significant role in the large 
scale structure of the universe. Finally, the neutrino asymmetry in early universe may be 
responsible for leptogenesis and then baryogenesis, etcetera. 

In the flat space, the neutrino oscillation is due to difference in rest masses between two 
mass eigenstates. However, Gasperini [l] flrst pointed out that presence of gravitational 
fleld affects different neutrino flavors differently which violates equivalence principle and 
thus governs oscillation, even if neutrinos are massless or of degenerate mass. Subsequently, 
the effect of strong gravitational fleld on neutrino oscillation was discussed {2!. Later it was 
shown that the neutrino oscillation with LSND data [.3| can be explained by degenerate or 
massless neutrinos with flavor non-diagonal gravitational coupling. It was further argued ^ 
that in weak gravitational fleld the flavor oscillation is possible with the probability phase 
proportional to the gravitomagnetic fleld. The oscillation was also shown to be feasible when 
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the maximum velocities of different neutrino differ each other, even if they are massless . 

All the above results are for flavor oscillation or/and without any rigorous effects of 
general relativity. However, as we show below explicitly that the nature of curvature plays 
very important role to determine whether the oscillation occurs or not, even under strong 
gravity. Therefore, the goal of the present paper is to address the neutrino— antineutrino 
oscillation, which violates lepton number conservation, focusing on the nature of space-time 
curvature and its special effect. In this case, any oscillation depends on the gravitational 
coupling strength. The main message is something new and unique compared to that in 
earlier studies mentioned above and many other under curved space-time [6]. 

While the neutrino— antineutrino oscillation under gravity is an interesting issue on its 
own right, the present result is able to address two long-standing mysteries in astrophysics 
and cosmology as well as related nuclear physics: (1) Source of abnormally large neutron 
abundance to support the r-process nucleosynthesis in astrophysical site. (2) Possible origin 
of baryogenesis. 

It has recently been shown that under gravity, which may bring the CPT violating inter- 
action, energy level for the neutrino is split up from that of the antineutrino M, [S] . Similar 
interactions not preserving Lorentz and CPT symmetry were noticed earlier |9(] in different 
contexts. Now the splitting of energy level for neutrinos from that of corresponding anti- 
neutrinos may create a difference in their evolutionary phases. This CPT violation effect 
is responsible for the transition over neutrino mass eigenstates, if such states which are 
the linear combination of a neutrino and an antineutrino state are possible to construct, in 
presence of a mechanism to break the lepton number conservation. 

The possible violation of local Lorentz symmetry in the kaon sector of the standard model 



lOj. It was also pointed 



was studied with both CPT conserving and non- conserving cases 

out ll| that in a CPT non-conserved case the neutrino and the antineutrino acquire different 



mass and thus there is a possibility of oscillation between them. Later, it was shown 12|] in a 
CPT non-conserved case that in absence of neutrino decay a neutrino evolves differently from 
an anti-neutrino of the same flavor and thus the oscillation between them. Recently, it has 
been argued |l3] that a linear superposition of two opposite helicity states can be described 
as an eigenstate of the de Sitter Casimirs which actually gives rise to the Majorana state. 
However, all these works neither explain the explicit physical origin of CPT violation nor 
demonstrate the dependence of oscillation upon the physically observable quantities. In the 
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present paper, we consider realistic situations when neutrinos are propagating in a strong 
gravitational field, e.g. the space-time close to a compact object or in the early universe 
era, such that the spin of a neutrino couples with the spin-connection to the space-time and 
thus generates a gravitational interaction even if there is no other interaction present. 

In the next section, we describe the basic formalism of the problem starting from the 
fermion Lagrangian density and then how does this generate the oscillation between a neu- 
trino and an antineutrino state. Subsequently, we discuss about the oscillation phase and the 
oscillation length in various natural space-times in §§III,IV. In §V, we discuss experimental 
bound on CPT violation. Finally we summarize the results in §VI. 



II. FORMALISM 



A. Fermion Lagrangian density in curved space-time 

We consider a locally flat coordinate system where the gravitational interaction comes in 
as an effective interaction. Then the neutrino Lagrangian density under gravity is 

m 

C = v^? [{^r^a -m)+ n'B,] ^ = Cf + Ci, (1) 

where 

5' = e'^^'^'^e,, (9,e^ + r^^eX) , e^, = g"^ , (2) 

where the choice of unit is, c = K = ks = 1- Naturally, the Lagrangian density in a local 
flat space-time can be split out into two parts. One is the free part similar to the flat space 
Lagrangian density apart from a multiplicative factor ^/g, and other is the gravitational 
interaction part which comes in due to an effective extension of a flat-space description. 



B. CPT status of the Lagrangian 

Ci may be a CPT and Lorentz violating interaction if the background curvature coupling, 
Ba, is constant in the local frame or CPT even [s^. Earlier, similar interaction terms were 



also considered in CP' 



origin in nature 



violating theories and string theory without any explicit possible 
If Ba does not flip in sign under CPT transformation, then Cj 



violates CPT. Actually, as described in detail in our previous work [15|, under the CPT 



transformation, associated axial-vector or pseudo-vector ('ip'j'^'j^ilj) changes sign. If Ba is 
treated as a background field in a local frame, then the interaction violates CPT explicitly. 
However, in the present case, with its functional form we can determine the explicit CPT 
status of Ba itself at a space-time point {t,x,y,z). If Ba{—x, —y, —z, —t) = —Ba{x,y,z,t), 
then Ba is CPT odd. On the other hand, CPT even Ba corresponds to Ba{—x, —y, —z, —t) = 
Ba{x,y, z,t). If Ba{x,y, z,t) is not an odd function under CPT [Ba{—x,—y,—z,—t) ^ 
—Ba{x,y,z,t)], then Cj comes out to be a CPT violating interaction along the space-time. 
It is the nature of background metric which determines whether Ba is odd or even under 
CPT and then the overall CPT status of the interaction ^. 

It is very important to note, and can be easily verified, that the gravitational coupling, 
Ba, is zero when the space-time is spherically symmetric and the metric does not consist of 
any off-diagonal element, e.g. for the Robert son- Walker universe, the Schwarzschild black 
hole, etcetera. Nevertheless, there may exist the spin-orbit coupling when the helicity of a 
neutrino is different from its chirality. However, the spin-orbit coupling effect is expected to 
be small compared to the effect due to Ba arised in a non-spherical symmetric space-time. 

C. Energy dispersion relations for neutrinos and anti-neutrinos 

Now we stick to the standard model so that a neutrino is solely left-handed and an 
antineutrino is solely right-handed. Therefore, when the background gravitational field is 
constant or Ba is CPT even, Cj, in Weyl's representation, can be written explicitly in 
combination of the neutrino and the antineutrino fields as 



^l^r7'^Ba = [ipnri^R ' i^L^^^L) Ba- (3) 
Then, the dispersion energy relations for the neutrino and the anti-neutrino are 



E, = ^J{p-By + m^ + Bo, ETj=^{p + By + m^-Bo. (4) 

As El, 7^ E-jj ^, the time evolution of a neutrino is different compared to that of an antineu- 
trino. 



^ The associated axial-vector, tp'^^-'^^ip, always changes sign under CPT. 

^ This is clearly due to presence of background curvature which appears as gravitational four vector Ba- 
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Let us now consider the Majorana neutrino when the anti-particle is basically the charged 
conjugated particle, and thus one can construct the Majorana spinor which has only left- 
handed component such that ip = ipi and ip'^ = ipl = ipji, which gives rise to the lepton 
number violating mass terms in C Therefore, eqn. ([T]) can be re-written explicitly for 
Majorana neutrinos as 



-9 



Thus, the Lagrangian may break the lepton number conservation and CPT both simultane- 
ously, depending of the background, which may lead to oscillation. 

D. Oscillation between neutrino and antineutrino states 



Now motivated by the neutral kaon system, we consider two distinct orthonormal eigen- 
states \Ey > and \Eu >, for a neutrino and an antineutrino type respectively, of the same 
flavor. Further we introduce a set of neutrino mass eigenstates at t = as 



|mi >= cos6 \Ei, > +sin9 \Etj >, |m2 >= —sinO \Ey > +cos9 \E^ > 



(6) 



Therefore, the oscillation probability for \mi{t) > at t = to \m2{t) > at a later time t = ti 
can be found as 
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[-sin6 <E^\+ cose < Etj\ ] [e'"'^''^^ cosO \E^ > +e-'^^'sin9 \Etj > 
sin^20 



(7) 



where for ultra-relativistic neutrinos, 

{E,y — Ejj)ti 



[Bo-m 



(8) 



The term within parenthesis in the right hand side of eqn. ([8]) is arised due to gravity, while 
the other term is the fiat space contribution which is zero when rest masses of the neutrino 
and the antineutrino are same and/or in the massless limit. In rest of the description we 
set Am^ = 0, as our aim is to establish any gravity effect which is expected to contribute 
dominant ly. 

Therefore, a non-zero oscillation phase, due to Ba 7^ 0, between a neutrino and an anti- 
neutrino is expected to appear which is gravitational in nature. The interesting point to 



note is that the present mechanism explicitly affects the neutrino— antineutrino pair unlike 
the cases considered earlier (e.g. 1^) when the gravity would mainly affect different 
ffavors. Note that a general field theoretic description of neutrino oscillation including the 
possible neutrino-antineutrino oscillation has already been worked out and has been shown 
that CPT violation provides oscillations without the additional sterile neutrinos [16;]. In 
a similar fashion, our mechanism also explains neutrino-antineutrino oscillation without 
incorporating sterile neutrinos. In general, it is the non-zero coefficient for CPT violation 
in the effective Lagrangian which generates neutrino-antineutrino oscillation in both the 
models, while we explicitly show that the CPT violation coefficient may originate from the 
background curvature. Both the models mainly emphasize the possibility of oscillation in 
absence of neutrino mass difference. The unconventional energy dependence in oscillation 
phase^ (for massless neutrinos) given by eqn. ([8]) also reflects its similarity to the earlier 



work 



16|. 



III. OSCILLATION PHASE IN SPECIFIC SPACE-TIMES 

It is very clear from eqns. ([1]), (jl]) and ([5]) that due to gravity the neutrino and the 
anti-neutrino acquire different effective mass which finally leads to a non-zero oscillation 
probability in two mass eigenstates. Therefore, in curved backgrounds, e.g. in the anisotropic 
phase of early universe when the space-time is non-ffat, close to a rotating compact object, 
the neutrino-antineutrino oscillation exists in presence of a suitable mechanism to violate 
lepton number conservation. Below we describe the oscillation according to specific space- 
times. 



A. Early Universe 

First we consider the anisotropic phase of axially symmetric early universe, when the 
GUT processes lead to lepton number non-conservation. With a simplified version of the 
Bianchi II model, 

ds^ = -de + S{tf dx^ + R{tf [dy^ + f{yf dz^] - S{tf h{y) [2dx - h{y) dz] dz, (9) 
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where f{y) = y and h{y) = —y'^/2. The corresponding orthogonal set of non-vanishing 
tetrad (vierbien) components can be chosen as 

e? = 1, el = f{y)R{t)S{t)/^f{yyR{ty + S{tmyy, = R{t), 

^JfiyyUity + S{tfh{y)\ el = -S{tfh{y)/^f{yYR{tY + S{tYh{yY. (10) 



el 



d 



We then obtain the components of B 

AR^S + ^y^RS^-2yS' 

2 {AyR'-^RS-y^S'){RS' -RS) , 
^= 8R^ + 2ym-S- (^^^ 

It is very clear from above that B^ and B^ do not flip sign under space-inversion, i.e. for 
y — > —y. Thus, it is not an odd function over the space-time for any of the Bianchi models 
and the form of Ba is such that Ba{—x, —y, —z, —t) ^ ±Ba{x, y, z, t). An example of a space- 



time where Bo{—x,—y,—z,—t) = —BQ{x,y,z,t) can be found in [l7j such that £/ under 
CPT does not change sign overall. Therefore, according to the discussion in §11. B, Ba leads 
to CPT violation at any space-time point {x,y,z,t). Hence, one can obtain the oscillation 
phase 6 from eqn. ([H]). If we consider a special case when S{t) = arbitrary constant = k, 
then B^, B"^ —>■ ai t oo ^, which verifies that the curvature coupling is insignificant in 
the present universe to exhibit any such oscillation. It also can easily be checked that in 
a spherically symmetric case, say for Robertson- Walker universe, B'^ vanishes at all t and 
thus the oscillation probability is zero. 

B. Around black holes 

Second example can be given for a space-time around the compact object when the 
curvature effect is significant. In early, the flavor oscillations around rotating black holes (in 

jj]. The black hole space-time, namely the Kerr 

2ffl as 



case of active galactic nuclei) were studied 

geometry, can be described in Cartesian-typed coordinate (t, x, y, z 

-2a 



ds = rjij dx'^ dx^ — 



Si Vj + O? Vi Vj 

p 



dx' dx^ (12) 



It is easy to compute that R{t) ~ (At — B)^^'^, when A and B are arbitrary constants [l8|. 
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where 



V2Mr 



2^2 



2 2 I 

a = , p =r + 

p 



a z 



ay -ax \ ( rx ry z\Jr'^ + a^\ 



gj2 _^ ^2 ' 0^2 _^ ^2 ' y' \^ ' ^^2 _^ q2 ' ^7^2 _^ ^2 ' 

Here a and M are respectively specific angular momentum and mass of the Kerr black hole 
and r is positive definite satisfying, 



- (x^ + y'^ + z^ - a^) - a^z^ = 0. (14) 



The corresponding non-vanishing components of tetrad (vierbien) are given as 201] 
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S3V2, 
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el = l-yS3V3. (15) 

From eqn. ([2]) we obtain the gravitational scalar potential 

= eiA (^36^ - ^26^) + e2A (die', - ^36^) + e3A (^26^ - d^e^) = (16) 



pV2r3 



2 - z^ 



where p^ = 2r^ + — — y 
Similarly, one can obtain B^,B^,B^. Clearly, from eqn. ( fT6l) . B^{—t,—x,—y,—z) = 
B^{t,x,y, z); above B^ is a CPT even function which makes £7 a CPT odd interaction. 
Therefore, the neutrino-antineutrino oscillation is feasible and the oscillation phase, 6, can 
easily be computed. It is then a trivial job to check that the oscillation phase for a non- 
rotating (Schwarzschild) compact object is zero, as there B'^ itself vanishes, by setting a = 
in the expression of B'^. 



IV. OSCILLATION LENGTH AND ITS CONSEQUENCE 

The above discussion establishes the fact that the gravity induced neutrino-antineutrino 
oscillation is expected to occur provided the metric consists of one off-diagonal term at least. 
The amplitude of oscillation is zero at ^ = 0, 7r/2 and maximum at = 7r/4. From eqns. ([7]) 
and ([8]) the oscillation length. Lose, by appropriately setting dimensions, is obtained as 



where B = Bq — \B\ is expressed in GeV unit and the neutrino is moving in the speed of 
hght. 

A. Around black holes 

One of the most fruitful situations for an oscillation of this kind to take place is the inner 
region of an accretion disk around a rotating compact object. Therefore, for a compact 
object of mass M = Mg Mq, the oscillation length at an orbit, r p = x M and z = H M, 
of the disk, assuming _B ~ i?o, is computed from eqns. (fT6ll and (fTTl) 



This resn 
(NDAF) 



t has a very important implication to the neutrino dominated accretion flow 



22| |. Although the NDAF and the related supernova are described in the 
existing literature ignoring any gravity induced oscillation effects, it is expected that the 
oscillation at an inner edge of the disk to be influenced by gravity that affects the neu- 
trino sphere and then the accretion dynamics and outflow which are directly related to the 
corresponding prediction of supernova explosion. It appears from eqn. (ITS!) that at the 
inner accretion disk, when x < 10, Lose varies from a few factors to several hundreds of 
Schwarzschild radii, depending on the location and the thickness of the disk, for a fast spin- 
ning compact object. This is interesting as the accretion disk can be extended upto several 
thousands of Schwarzschild radii. 

One of its important consequences is to the r-process nucleosynthesis. Supernova is 
thought to be the astrophysical site of the r-process nucleosynthesis. During supernova, 
neutron capture processes for radioactive elements take place in presence of abnormally 
large neutron flux. However, how does the large neutron flux arise is still an open question. 
There are two related reactions: 

n + z/g^p + e", p + h'e^n + e'^. (19) 

If z/e is over abundant than u^, then from eqn. f|T9|) neutron production is expected to be 
more than proton production into the system. Therefore, the possible conversion of Uf, to 
due to the gravity induced oscillation explains the overabundance of neutron. 
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B. Early Universe 



If the oscillation is considered in early universe with anisotropic phase, then at the GUT 
scale B ~ 10^ GeV [15:]. From eqn. f[T7|) . this leads to Lose ~ 10~^^km which is 10^^ orders 
of magnitude larger than the Planck length. This has an important implication as the size of 
universe, x = 

Jocdt'/R{t'), at the GUT era is within ~ 10^^ times of the Planck. Therefore, 
the oscillation may lead to leptogenesis and then to baryogenesis by electro-weak sphaleron 
processes due to B — L conservation, what we see today. 



C. Atmosphere of Earth 

If we consider a case of neutrino— antineutrino oscillation in the atmosphere of earth 
where the curvature effect, B ~ 10"^'' GeV, can be found on a satellite orbiting earth with 
velocity f ~ 1 km/sec [17|, then Lose at that orbit comes out to be of the order of 10^*^ km. 
This is quite large that the satellite revolves about 10^^ times in the time taken to complete 
only one oscillation. This is because the neutrino— antineutrino oscillation is difficult to 
observe in the atmosphere of earth. 



V. BOUNDS ON CPT VIOLATION FROM EXPERIMENTS AND OUR MODEL 

Kostelecky and his collaborators 9| already argued that the CPT violation, if present, 
would be a very small effect and then they argued for very stringent constraints to the 
corresponding coefficients. The experiment E773 at Fermilab published a bound to mass 
difference between particle and antiparticle in the neutral kaon sector as = l^^o ~ 

n 

rrij^ol/rnxo < 1.3 x lO""*^^ [23]. It was also argued that any observable effects must be 



suppressed due to involvement of, perhaps, higher-level fields of Planck scale mass Mpi [24 1. 
When the electroweak scale mass niew ~ lO^GeV, rriew/Mpi ^ 10^^^ provides the natural 
dimensionless quantity which governs the suppression. In the neutral kaon system, one 
could expect ~ rriew/Mpi which is just below the present bound. The same authors also 
pointed out that rriew/Mpi might be expected to arise in a gravitational mechanism that 
violates CPT. 

In the present case, we have constructed the neutrino mass eigenstates under gravity 
described in §11. D in the spirit of the neutral kaon mass matrix when the mass difference 
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between neutrino and antineutrino is ~ 2B. Now from §IV, for black holes of mass range, 
e.g. IOMq < M, < lO^M©, the oscillation length comes out to be 10km ^ Lose ~ lO^km, 
(which is equivalent to a few factors to a few hundreds Schwarzschild radii, as mentioned 
in §IV.A), provided lO^^^GeV ^ B ^ lO^^'^GeV. Therefore, the necessary B to govern any 
significant oscillation in an accretion disk is very small and is quite compatible with the 
experimental bound. 

In the GUT scale, as described in §IV.B, the physically interesting oscillation length may 
range e.g. as lO'^'^km ^ Lose ~ lO^^km, provided lO^GeV ^ B ^ 10^^°GeV, when size of 
the universe is ^ lO^^^km. If we compare this B with the mass scale (mcuT ~ lO^^GeV) 
of the system, then it comes out to be one part in 10^^ to one part in 10^^. 

Much more stringent bounds can be put on the gravitational pseudo four- vector potential 
Ba (which is equivalent to B) from the tests for CPT and Lorentz violation [25]. Cj in non- 
relativistic limit is equivalent to s.B due to the interaction between the fermion spin s 
and the external field B. One can measure this interaction energy in experiments where a 
macroscopic number of fermions can be polarized in the same direction. For example, in the 
Eot-Wash II experiment 261], the spin-polarized torsion balance has = 8 x 10^^ aligned 
spins (with very negligible net magnetic moment). There will be a torque r = {N/'n-)AE 
on such a torsion balance where AE = \B\ is the difference in energy between the fermion 
spins polarized parallel and antiparallel to the external B field. From this experiment 26| . 
one can measure up to |5| ~ 10~^'^GeV {2?]. The field B also can be probed by measuring 
the net magnetization in a paramagnetic material using a squid 28|]. An external field B 
appears as an effective magnetic field of strength i?efT = {B/ ^b)- The magnitude of B^q can 
be probed in this experiment as |-Beff| = lO^^^G which reflects B ~ lO^^^GeV 27, [2^. 



VI. SUMMARY 



Although the flavor neutrino oscillation has already been described in several occasions 
successfully, the underlying mechanism for the neutrino— antineutrino oscillation is still not 
well understood. While in the first case, lepton number of the system remains conserved, 
the second case is a lepton number violating phenomenon that causes it more difficult to 
explain and establish. We have argued that in presence of space-time curvature violating 
CPT symmetry, the energy level of a neutrino splits up from that of an antineutrino which 
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gives rise to their different effective mass. This affects their distribution and results their 
difference in phases of evolution. Then, motivated by the neutral kaon system, we have 
described two mass eigenstates as linear combination of a neutrino typed and an antineutrino 
typed state. Therefore, the oscillation between mass eigenstates is expected. 

We have explicitly demonstrated this lepton number violating oscillation with some nat- 
ural examples. We have shown that in the anisotropic phase of early Universe, when the 
space-time is non-flat, the neutrino— antineutrino oscillation is feasible which may lead to 
the leptogenesis and then to the baryogenesis in presence of a suitable mechanism to violate 
the lepton number conservation. This may be the possible origin of the matter-antimatter 
asymmetry in Universe what we see today. 

On the other hand, in the inner region of an NDAF where the general relativistic effect 
is important, this oscillation may affect the accretion dynamics and outflow. It has recently 



been argued [22| that the energy carried out in the outflowing wind from such disks could 
be as high as 10^^ erg which might be sufficient to convert a failed supernova explosion 
into a successful one. They have also noted that the neutron to proton ratio is large in 
several regions of the disk. The neutron-rich regions are thought to be the astrophysical 
site for r-process nucleosynthesis. However, it is unclear that how does the large neutron 
abundance over proton arise in the first place. Our calculation suggests that this is due to 
larger antineutrino, due to oscillation, over the neutrino which favors the proton to decay 
into neutron compared to the neutron to proton. One now needs to study the NDAF in detail 
incorporating oscillation effects and to see how does the disk structure including neutrino 
sphere change and then affect the energy budget. 
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